function [Z P K] = ZPKadjust(z, p, G)
%G es la ganancia deseada

    %p, z como filas
    dims = size(z);
    if dims(1) > dims(2)
        z = z.';
    end

    dims = size(p);
    if dims(1) > dims(2)
        p = p.';
    end

    n = length(p);
    P = p(imag(p)>1e-6);
    pconj = conj(P);
    if ~rem(n,2) % Si 'n' es par
        P = [P pconj];
    else
        P = [P pconj real(p(round(imag(p)*1000)/1000 == 0))];
    end  
    
    n = length(z);
    if(n ~= 0)
        Z = z(imag(z)>0);
        zconj = conj(Z);
        if ~rem(n,2) % Si 'n' es par
            Z = [Z zconj];
        else
           Z = [z(imag(z) == 0); Z; zconj];
        end
    else Z = [];    
    end
    
    
    K = 10^(G/20) * abs(prod(p)) / abs(prod(z));

end
